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Expression Analysis in Multiple Muscle Groups and Serum 
Reveals Complexity in the MicroRNATranscriptome of 
the mdx Mouse with Implications for Therapy 

Thomas C Roberts 1 , K Emelie M Blomberg 2 , Graham McClorey 1 , Samir EL Andaloussi 1 2 , Caroline Godfrey 1 , Corinne Betts 1 , 
Thibault Coursindel 3 , Michael J Gait 3 , CI Edvard Smith 2 and Matthew JA Wood 1 



MicroRNAs (miRNAs) are a class of small RNAs that regulate gene expression and are implicated in wide-ranging cellular 
processes and pathological conditions including Duchenne muscular dystrophy (DMD). We have compared differential miRNA 
expression in proximal and distal limb muscles, diaphragm, heart and serum in the mdx mouse relative to wild-type controls. 
Global transcriptome analysis revealed muscle-specific patterns of differential miRNA expression as well as a number of 
changes common between tissues, including previously identified dystromirs. In the case of miR-31 and miR-34c, upregulation 
of primary-miRNA transcripts, precursor hairpins and all mature miRNAs derived from the same transcript or miRNA cluster, 
strongly suggests transcriptional regulation of these miRNAs. The most striking differences in differential miRNA expression 
were between muscle tissue and serum. Specifically, miR-1, miR-133a, and miR-206 were highly abundant in mdx serum 
but downregulated or modestly upregulated in muscle, suggesting that these miRNAs are promising disease biomarkers. 
Indeed, the relative serum levels of these miRNAs were normalized in response to peptide-phosphorodiamidate morpholino 
oligonucleotide (PMO) mediated dystrophin restoration therapy. This study has revealed further complexity in the miRNA 
transcriptome of the mdx mouse, an understanding of which will be valuable in the development of novel therapeutics and for 
monitoring their efficacy. 
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INTRODUCTION 

Duchenne muscular dystrophy (DMD) is an X-linked, reces- 
sive, fatal disorder caused by mutations in the gene that 
encodes the dystrophin protein. 12 DMD is the most common 
muscular dystrophy and affects 1 in 3,500 live male births. 
Disease progression is characterized by progressive myofi- 
ber degeneration, muscle wasting, fibrosis and loss of ambu- 
lation. Degeneration of cardiac muscle and the diaphragm 
results in cardiac and/or respiratory failure with subsequent 
fatality in the second or third decade of life. 

Although the primary genetic defect in DMD has been 
known for several decades, effective treatment has been lim- 
ited by a lack of understanding of the complex, multifactorial 
disease pathophysiology. Loss of functional dystrophin at the 
sarcolemma and concomitant disruption of the dystrophin- 
associated protein complex leads to wide-ranging cellular 
pathologies. These include: myofiber necrosis, impaired sig- 
naling activities (e.g., NO signaling via neuronal nitric oxide 
synthase), 3 4 loss of sarcolemmal integrity, 5-7 irregular calcium 
homeostasis, 5 8 " 10 increased oxidative stress, 11 chronic cycles 
of degeneration and regeneration leading to satellite cell 
depletion, 12 and inflammation. 13 The utility of gene expres- 
sion profiling has the potential to identify the sequence of 



events that leads from dystrophin loss to these pathological 
processes. 

MicroRNAs (miRNAs) are a class of small (21-25 nt) 
non-protein-coding RNAs that regulate gene expression 
by a homology dependent mechanism distinct from, but 
related to the small interfering RNA-mediated RNA interfer- 
ence pathway. 14 Primarily, miRNAs function by binding to the 
3' untranslated regions of target mRNA transcripts forming 
imperfect RNA duplexes whereby they facilitate either trans- 
lational repression or decapping and subsequent endonu- 
cleolytic degradation of the transcript. 15-17 A 3' untranslated 
region may contain many miRNA target sites and a single 
miRNA can target multiple mRNA transcripts. The majority 
of mammalian mRNAs (-60%) contain conserved miRNA 
targets underlining the importance of miRNA regulation. 18 
Consequently, miRNAs act to modulate global gene expres- 
sion through the regulation of families of related transcripts. 19 
Many miRNAs show tissue-specific expression 20 and a 
number of miRNAs are primarily expressed in skeletal and 
cardiac muscle (e.g., miR-1, miR-1 33, and miR-206). 2122 
Transcription of muscle-specific miRNAs has been shown 
to be regulated by transcription factors involved in muscle 
growth and development (e.g., MYOD, MEF2, and MYF5). 23-25 
In addition, non muscle-specific miRNAs (miR-26a, miR-1 81, 
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and miR-486) have also been found to regulate skeletal mus- 
cle processes. 26-28 

The discovery that unique miRNA expression patterns were 
associated with different tumour types 29 has led to investiga- 
tion of miRNA profiles in a wide variety of other diseases. 
Two studies profiling DMD-associated miRNAs have identi- 
fied common differentially expressed miRNAs despite con- 
siderable differences in results. 3031 (Supplementary Table 
S1). For example, Greco et al. showed a ~70-fold increase in 
miR-31 expression in both mdxan6 human muscle whereas 
Eisenberg et al. did not identify this miRNA as being differ- 
entially expressed (in DMD patient biopsies). Conversely, the 
most differentially expressed miRNA in the Eisenberg study, 
miR-146b, was not significantly changed in the study by 
Greco et a/. The differences in results between these studies 
may be due to the methodologies used (microarray and Taq- 
Man array, respectively) and the assaying of different muscle 
groups. 

Here, we have investigated miRNA expression in a range 
of muscle groups and serum in mdxm\ce compared with age- 
and sex-matched wild-type controls. We have used miRNA 
microarray technology to perform an initial comparative 
miRNA analysis between quadriceps femoris, diaphragm, 
and heart muscles. Using small RNA TaqMan assays we 
subsequently validated a number of differentially expressed 
miRNAs identified on the arrays and expanded the study to 
include the triceps and tibialis anterior (TA) muscles, and 
serum. We have focused on a set of 11 miRNAs which we 
refer to as "dystromirs."This term has been used previously 32 
but here we define it more specifically as miRNAs which 
have been previously found to be differentially expressed in 
DMD patients and/or dystrophic animal models. The dystro- 
mirs relevant to this study are described in Table 1 . 

A comprehensive understanding of miRNA expression 
in the mdx mouse is important because: (i) miRNAs are 
promising biomarkers for noninvasive monitoring of disease 
progression and response to experimental therapies, (ii) it 
contributes to the understanding of the series of molecular 
events leading from loss of dystrophin to dystrophic pathol- 
ogy, and (iii) it forms the basis for the development of novel 
miRNA-based therapeutics. The present study is the first to 
compare differential miRNA expression in the mdx heart and 
diaphragm and has directly compared multiple muscles and 
serum from the same set of animals. The findings point to 
an unexpected layer of complexity in the mdx mouse miRNA 
transcriptome. 



RESULTS 

miRNAs microarray analysis 

To identify differentially expressed mature miRNAs associ- 
ated with the DMD phenotype, miRNA microarray analy- 
sis was performed on quadriceps, diaphragm and cardiac 
muscles harvested from 8-week-old C57BI/10 (wild-type) 
and mdx mice. Statistically significant (P < 0.05) changes in 
mature miRNA expression were identified across all of the 
tissues; 188 miRNAs in quadriceps (false discovery rate 
(FDR) = 6%), 181 miRNAs in diaphragm (FDR = 5%) and 
64 miRNAs in heart (FDR = 26%). Analysis of the full miRNA 
transcriptome reveals widespread differences between 



Table 1 Dystromirs of interest in this study 



miRNA 



Function 



Evidence 



miR-1 



Myogenesis 



miR-1 33a Myogenesis 



miR-206 Myogenesis 



miR-21 Fibrosis 
miR-29c Fibrosis 



miR-31 Regeneration 



Regulates dys- 
trophin 



miR-34c 



Regeneration/ 
degeneration 



miR-1 46b Myogenesis 



miR-1 99a-3p Unknown 
miR-221 Myogenesis 



miR-223 Inflammation 



Downregulated in mdx adductor 
and DMD patients 30 

Downregulated in mdx TA and 
CXMDj TA 39 

Downregulated in mdx 
gastrocnemius 43 

Highly abundant in mdx serum and 
DMD patient serum 32 - 33 

miR-1 promotes muscle 
differentiation 21 

Downregulated in mdx TA and 
CXMDj TA 39 

Downregulated in mdx 
gastrocnemius 43 

Highly abundant in mdx serum and 
DMD patient serum 3233 

miR-1 33a promotes muscle 
proliferation 21 

Upregulated in mdx adduc- 
tor and diaphragm and in DMD 

patients 303137 

Upregulated in mdx TA and 
downregulated in CXMD d TA 39 

Upregulated in mdx gastrocne- 
mius 43 

Highly abundant in mdx serum and 
DMD patient serum 3233 

miR-206 promotes muscle 
differentiation 22 

Upregulated in DMD patients 31 . 

miR-21 promotes fibrosis 50 

Downregulated in mdx adductor 
and in DMD patients 3031 

Downregulated in mdx 
gastrocnemius 43 

miR-29 regulates collagens and 
elastin 51 

Expression of miR-29 in mdx 
gastrocnemius reduces fibrosis 43 

Upregulated in mdx adductor and 
DMD patients 30 (most upregulated 
miRNA in this study). 
miR-31 is strongly induced in 
ischemia damaged myofibers 30 

miR-31 regulates dystrophin 
expression 38 

Upregulated in mdx adductor and 
DMD patients 30 

miR-34c is strongly induced in 
ischemia damaged myofibers 30 

miR-34c promotes cell cycle 
withdrawal and apoptosis 52 

Upregulated in DMD patients 31 
(most upregulated miRNA in this 
study). 

miR-1 46 promotes proliferation 
of C2C12 myoblasts 53 
Upregulated in DMD patients 31 
Upregulated in DMD patients 31 . 
miR-221/222 is highly upregu- 
lated in terminally differentiated 
myoblasts 54 

Upregulated in mdx adductor 
and DMD patients 30 
Upregulated in mdx 
gastrocnemius 43 
Possible immune infiltrate 55 



Abbreviations: DMD, Duchenne muscular dystrophy; miRNA, microRNA; 
TA, tibialis anterior. 
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these tissues (Figure 1a). More miRNAs were differentially 
expressed in quadriceps and diaphragm compared to heart 
and, in addition, the observed changes in the heart were 
generally of lesser magnitude than in quadriceps and dia- 
phragm (Figure 1b). Overlap between miRNAs commonly 
differentially expressed across all three tissues is indicated 
in the Venn diagram (Figure 1c). For each tissue there was 
a set of miRNAs found to be differentially expressed only in 
that tissue suggesting tissue-specific differences in patho- 
logical processes. 1 8 miRNAs that were significantly differen- 
tially expressed in all three tissues are shown graphically in 
Figure 1d and fold changes listed in Supplementary Table 
S2. In addition to the dystromirs listed in Table 1, miR-1944, 
miR-296-3p, and miR-543 showed increased expression, 
and miR-185, miR-151-5p, miR-151-3p, miR-345-3p, miR- 
139-5p, miR-143, and miR-378* showed decreased expres- 
sion in all three mdx tissues. To our knowledge these miRNAs 
listed above have not been previously identified as differen- 
tially expressed in miRNA profiling studies of the mdx mouse 
or dystrophic patients. 

RT-qPCR validation of array data 

In order to validate the results of the miRNA microarray, 
reverse transcription-quantitative PCR (RT-qPCR) analy- 
sis was performed on the same RNA samples used for the 
arrays. Small RNA TaqMan assays designed to detect only 
the mature miRNA species were used for all of the dystro- 
mirs in Table 1. In order to determine whether inferences 
drawn from array data can be applied to other tissues, the 
study was expanded to include dystromir expression in addi- 
tional tissues; triceps and TA muscles (Figure 2a-k and fold 
changes listed in Supplementary Table S3). Serum was 
also included as changes in serum miRNA abundance have 
recently been reported in the mdx mouse, CXMDj dog and 
DMD patients. 3233 miR-1 and miR-133a were downregulated 
or unchanged in all mdx muscles whereas in serum both 
miRNAs were increased by 43.6-fold and 50.6-fold, respec- 
tively (Figure 2a,f). miR-206 was increased in all mdx groups 
measured with a 52.9-fold increase observed in serum (Fig- 
ure 2i). miRNA-21, -31, -34c, and -221 were increased in 
all muscles except the heart (only modest increases were 
observed in mdx serum in each case) (Figure 2b,d,e,j). miR- 
223 was found to be increased in the skeletal muscles but 
not in the diaphragm, heart, or serum (Figure 2k). Notably, 
miR-1 46b was only significantly increased in mdx quadriceps 
(Figure 2g) and similarly, miR-1 99a-3p was only increased in 
mdx quadriceps and diaphragm (Figure 2h). 

Generally, the RT-qPCR validation closely matches the 
miRNA microarray data (Figure 3a,b and Supplementary 
Table S4). Hierarchical clustering analysis was used to group 
fold changes for the respective miRNAs and tissues. As 
expected, quadriceps, TA, and triceps muscles cluster away 
from diaphragm and heart. Similarly, all the tissues cluster 
away from serum (Figure 3b). In quadriceps and diaphragm 
fold changes were highly similar between the two methodolo- 
gies used (e.g., in the case of miR-31 expression; 70.9-fold by 
RT-qPCR and 74.3 by microarray) with few exceptions. A num- 
ber of miRNAs were changed in the same direction between 
the two methodologies but did not reach significance at the 
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Figure 1 MicroRNA (miRNA) transcriptome analysis. Quadri- 
ceps, diaphragm, and heart tissue from 8-week-old C57/BI10 [wild 
type (WT)] and mdx mice {n = 4) were harvested and miRNA mi- 
croarray analysis performed, (a) Heatmap of whole miRNA tran- 
scriptome of all three tissues for all miRNAs with a statistically 
significant (P< 0.05) change in expression between WT and mdx 
mice, (b) Scatter plot showing the relative spread of fold changes 
across all three tissues. The most upregulated and most down- 
regulated miRNAs across all of the arrays are highlighted, (c) Venn 
diagram showing common changes in miRNA expression between 
the tissues, (d) Heatmap showing miRNAs with statistically signifi- 
cant expression changes in all three tissues. Hierarchical clustering 
analysis was used to identify patterns of miRNA expression in the 
mdx quadriceps, diaphragm, and heart. Red signifies an increase 
in expression and green, a decrease in expression. Heatmap scale 
bars indicate log 2 -fold changes. 

95% level as measured by RT-qPCR (i.e., miR-1, miR-29c, 
and miR-1 33a in quadriceps, and miR-29c, miR-1 46b, miR- 
199a-3p, and miR-223 in diaphragm). Interestingly, miR-1 
was unchanged in the diaphragm array dataset but was sig- 
nificantly decreased as measured by RT-qPCR (-2.27-fold). 
In the heart dataset only miR-206 was significantly increased 
in mdx heart as measured by RT-qPCR but not by microarray. 
No other miRNAs were significantly changed in mdx heart as 
measured by RT-qPCR. 
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The expression of one miRNA (miR-29c) did not match the 
microarray data. According to the miRNA microarray miR-29c 
was increased in mdx heart (1.29-fold, P= 0.0373) although 
the significance corresponds to a false discover rate of 0.24 
and the possibility of a false positive cannot be excluded. 
Similarly, miR-29c expression was not significantly different 
as measured by RT-qPCR (Figure 2c and Supplementary 
Table S4). 

In general, the cycle threshold (Ct) values for all miRNAs 
measured fell within the range of 20-30 cycles which was 



comparable to the reference miRNAs (Ct values around 
25-26). miR-31 and miR-34c were expressed at low levels 
(Ct values 30-35) in all the C57 tissues. Raw fluorescence 
values from the miRNA microarray also indicate low levels 
of expression for these miRNAs (data not shown). miR-206 
was present at very low levels in all of the heart samples (Ct 
values -35) which is unsurprising given that expression of 
miR-206 is generally restricted to skeletal muscle. 22 

Interestingly, the abundance of dystromirs in the mdx 
serum did not match the pattern of expression observed 
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in the tissues. This was reflected in the clustering analy- 
sis which showed that serum clustered separately from 
the tissues (Figure 3b). In particular, miR-1 and miR-133a 
were increased in mdx serum and were either decreased 
or unchanged in the mdx muscle tissues. Similarly, miR- 
206 was highly increased in serum but only moderately 
increased in the muscle tissues. Conversely, miR-31, was 
increased ~50-70-fold in quadriceps, TA and triceps but 
was only 2.3-fold increased in serum. Other dystromirs 
including miR-21, miR-34c, miR-221, and miR-223 followed 
similar patterns of differential expression between serum 
and muscle. 

Transcriptional dystromir upregulation 

Mature miRNAs are derived from hairpin structures 
residing within long primary transcripts. Often multiple 
miRNA hairpin structures are found in a single primary 
transcript, so-called miRNA clusters. After cleavage by 
Dicer, only one arm of the miRNA hairpin is loaded into 
the miRISC complex. miRISC loading is asymmetric and 
depends on the relative thermodynamic stability of the 5' 
ends of each strand (i.e., the strand with most unstable 
5' end is preferentially loaded). Inspection of the miRNA 
microarray data revealed that, in many cases, all of the 
miRNAs from a primary transcript were concordantly dif- 
ferentially expressed. For example, miR-31 and miR-31* 



(the minor miRNA species derived from the 3' arm of the 
miR-31 hairpin) were both increased in mdx quadriceps 
and diaphragm. Similar results were seen for the major and 
minor species of other dystromirs; miR-21, miR-1 46, and 
miR-1 99a (Supplementary Table S5) and for other non- 
dystromirs; miR-1 99b, miR-362, miR-675, miR-501, 
and miR-532 (data not shown). Similarly, the miRGen 
resource 34 was utilized to identify miRNA clusters within 
the mouse genome that are differentially expressed in the 
mdx mouse. Two differentially expressed dystromir-con- 
taining clusters were identified: (i) the miR-34 cluster (i.e., 
mature miR-34c, miR-34c*, miR-34b-5p, and miR-34b-3p 
all upregulated) and (ii) the miR-221/miR-222 cluster (i.e., 
mature miR-221 and miR-222 upregulated). A number of 
other differentially expressed miRNA clusters were identi- 
fied and are listed in Supplementary Table S6. Notably, 
mature miRNAs from the miR-1 88/miR-532 and miR-296/ 
miR-298 clusters were concordantly upregulated in both 
mdx quadriceps and diaphragm. 

The miRNA microarray chips also contained probe sets 
that detect pre-miRNA hairpins. Twelve pre-miRNA hair- 
pins were significantly (P < 0.05) differentially expressed 
across all three tissues (Figure 4a). Each pre-miRNA 
was differentially expressed in at least two tissues and 
two pre-miRNAs, pre-miR-181a-2 and pre-miR-185, were 
differentially expressed in all three tissues. Expression of 




Figure 2 Reverse transcription-quantitative PCR (RT-qPCR) analysis of dystromir expression. Tibialis anterior (TA), triceps, quad- 
riceps, diaphragm, heart, and serum from 8-week-old C57/BI10 [wild type (WT)] and mdx mice were harvested and small RNATaqMan 
reverse transcriptase-quantitative PCR (RT-qPCR) performed to determine relative expression of the following dystromirs: (a) miR-1, (b) 
miR-21, (c) miR-29c, (d) miR-31, (e) miR-34c, (f) miR-1 33a, (g) miR-1 46b, (h) miR-1 99a-3p, (i) miR-206, (j) miR-221, (k) miR-223. Gray 
bars represent WT samples and black bars represent mdx samples. Relative fold changes were determined by the Pfaffl method. Tissue 
microRNA (miRNA) expression was normalized to miR-1 6 and serum miRNA expression normalized to miR-223 (consequently serum 
expression of miR-223 is omitted from panel (k)). All values are mean + SEM. *P < 0.05, **P < 0.01, ***P < 0.001. For tissue samples 
n = 4. For WT serum, n = 3 and for mdx serum n = 4. 
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Figure 3 Heatmaps of mcfxdystromir expression. Heatmaps of dystromir expression as determined by (a) microRNA (miRNA) microar- 
ray analysis and (b) small RNA TaqMan RT-qPCR. Hierarchical clustering analysis was used to identify patterns of miRNA expression 
across the mdx tissues and serum. Red signifies an increase in expression and green, a decrease in expression. Heatmap scale bars 
indicate log 2 -fold changes. 
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Figure 4 Differential expression of precursor-miRNAs and primary-miRNAs. (a) Statistically significant changes in precursor-miRNA 
(pre-miRNA) expression in mdx quadriceps, diaphragm, and heart (n = 4) as determined by microarray analysis. Hierarchical clustering 
analysis was used to identify patterns of pre-miRNA expression in the mdx quadriceps, diaphragm, and heart. Red signifies an increase in 
expression and green, a decrease in expression. Heatmap scale bars indicate fold changes, (b) Changes in primary-miRNA expression in 
mdx quadriceps for pri-miR-31 , pri-miR-34c and pri-miR-206 as measured by RT-qPCR. Relative values were determined using the relative 
standard curve method and primary-miRNA expression normalized to PpiB. Values are mean + SEM, n = 4, * P < 0.05, ** P < 0.01 . 



both pre-miR-31 and pre-miR-34c was increased in mdx 
quadriceps and diaphragm which matches the observed 
increases in mature miR-31 and miR-34c species, 
respectively. 

Increases in mature miRNA levels could be due to either: 
(i) enhanced post-transcriptional Drosha/Dicer process- 
ing or (ii) upregulation of the primary-miRNA transcript. 
To discriminate between these possibilities we performed 



RT-qPCR with primer/probe assays that detect primary- 
miRNA transcripts for miRs-31, -34c, and 206. Statistically 
significant increases were detected for all three miRNAs in 
mdx quadriceps (Figure 4b). The increases in these pri- 
miRNAs correlated well with their mature miRNA counter- 
parts (correlation coefficient = 0.980). Taken together these 
data suggest that miR-31 , -34c, and -206 are upregulated at 
the level of transcription. 
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Figure 5 Pip6e-phosphorodiamidate morpholino oligonucleotide (PMO) treatment normalizes serum dystromir abundance. 

Twelve-week-old mdx mice were injected with a single 12.5 mg/kg dose of Pip6e-PMO intravenously and harvested 2 weeks later. Serum 
samples were analyzed for the expression of (a) miR-1, (b) miR-133a, and (c) miR-206 by small RNATaqMan RT-qPCR. miRNA levels 
were normalized to miR-223 expression and fold changes presented relative to the wild-type C57 average, (d) Representative immuno- 
fluorescence images of tibialis anterior muscle showing restoration of dystrophin at the sarcolemma in mdx mice following treatment with 
Pip6e-PMO. Samples were co-stained with laminin to indicate muscle fibers, (e) RT-PCR shows skipping of dystrophin exon 23 across 
various muscle groups. Values are mean ± SEM, n = 4, *P < 0.05. 



Pip6e-PMO treatment partially normalizes serum 
dystromir abundance 

Circulating miR-1, miR-1 33a, and miR-206 have been pro- 
posed as serum biomarkers for the progression of DMD. 32 - 33 
In order to determine whether the serum levels of these dys- 
tromirs responded to antisense oligonucleotide-mediated exon 
skipping therapy, 12-week-old mdx mice were treated with a 
single 12.5 mg/kg tail vein injection of Pip6e-conjugated phos- 
phorodiamidate morpholino oligonucleotide (PMO). (This con- 
jugate is a novel variant of Pip5e-PMO that had previously been 
demonstrated to induce skipping of exon 23 from the mature 
dystrophin transcript to restore restoration of dystrophin protein 
expression and muscle function 35 (also see C. Betts, A.R Saleh, 
A.A. Arzumanov, S.M. Hammond, C. Godfrey, T. Coursindel et 
al. related manuscript). Treated mice were harvested 2 weeks 
after injection and TA, quadriceps, heart and diaphragm mus- 
cles collected along with sera. Comparisons were made against 



age-matched wild-type and non-treated mdx controls. Treated 
animals showed a partial normalization of the high serum lev- 
els of miR-1 (18.1-3.6-fold), miR-1 33a (32.3-3.8-fold), and 
miR-206 (65.9-10-fold) observed in the mdx mouse (Figure 
5a-c). In all cases serum dystromir abundance was reduced 
to levels close to wild-type controls. Restoration of dystrophin 
protein expression by Pip6e-PMO treatment was confirmed in 
all muscles analysed by immunofluorescence (representative 
image of restoration in the TA shown in Figure 5d) and by the 
detection of mature dystrophin transcripts lacking exon 23, as 
determined by RT-PCR analysis (Figure 5e). 

DISCUSSION 

We have characterized differential miRNA expression 
between different muscle groups and serum of the mdx 
mouse. This study is the first to provide a complete profile by 
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microarray analysis of miRNA expression in the mdx heart 
and diaphragm, tissues which are central to the dystrophic 
pathology and disease progression in DMD patients. The 
use of miRNA microarray analysis has enabled the detec- 
tion of the maximum number of miRNAs possible (miRBase 
v15) whereas TaqMan Array Card technology is limited to 
381 miRNA assays. Generally, the mdx diaphragm and heart 
show less differential miRNA expression, both in number of 
differentially expressed miRNAs and in the magnitude of fold 
changes, reflecting a lower level of degeneration in the dia- 
phragm and heart in mice of this age. This is consistent with 
what is already known about the natural history of disease 
progression. 36 In addition, by assaying different tissues from 
the same animals we have shown that there are important 
differences between the skeletal muscle groups themselves. 
In general, quadriceps, TA and triceps muscles show similar 
levels of differential dystromir expression. However, significant 
changes in the expression of miR-146b and miR-199a-3p 
were not observed in all tissues. (Figure 2g,h and Supple- 
mentary Table S3). Differential expression between tissues is 
one explanation for contradictions between previous reports. 
For example, differential expression of miR-146b has been 
reported by some investigators 31 but not others 30 (Supple- 
mentary Table S1). One of the first studies of miRNA expres- 
sion (also in the 8-week-old mdx mouse) showed that miR-206 
expression was increased in the A77dxdiaphragm but not in the 
soleus or plantaris muscles. 37 The reported 4.5-fold increase 
in miR-206 expression in diaphragm is closely matched by 
our study (4.7-fold increase). However, the previous study 
suggested that miR-206 was upregulated in diaphragm but 
not in hindlimb on account of the diaphragm being the most 
severely affected mdx muscle. Our results do not support this 
conclusion as we have detected increased miR-206 expres- 
sion in two hindlimb mdx muscles (TA and quadriceps) and in 
the mdx forelimb (triceps). In all of these cases the increase 
in miR-206 expression is moderately higher than in the dia- 
phragm. To reconcile the differences between these studies 
we note that the study by McCarthy et al. utilized semiquantita- 
tive PCR and image quantification from PAGE gels, which is a 
much less quantitative and sensitive methodology than small 
RNA TaqMan RT-qPCR. Alternatively, a lack of increase in 
miR-206 expression in the mdx soleus and plantaris may rep- 
resent further complexity in the mdx muscle transcriptome. 

The most differentially expressed miRNA of all array data- 
sets was miR-31 (~70-fold increase). This is consistent with 
the study by Greco et al. which showed similar magnitude 
fold change in the adductor muscles of mice of the same age 
and in DMD patient biopsies (from the quadriceps femoris). 30 
Interestingly, Eisenberg et al. did not report miR-31 as being 
differentially expressed in a microarray study of patients with 
diverse muscle disorders (including DMD and Becker muscular 
dystrophy) 31 (Supplementary Table S1). We have observed 
high levels of miR-31 across all of the skeletal muscles stud- 
ied. This is encouraging as antagonism of miR-31 activity has 
recently been proposed as a synergistic therapeutic strategy 
to exon-skipping-mediated dystrophin restoration. 38 

Our results are consistent with those reported by Yuasa 
et al. who profiled miR-1, miR-133, and miR-206 in the TA 
muscle of mdxn\\ce and the CXMD d dystrophic dog model. 39 
Similarly, two recent studies have identified miR-486 as 



downregulated in mdx skeletal muscle. 2840 We have con- 
firmed this result as miR-486 expression was reduced 1 .3-fold 
(P = 0.0002) in the mdx quadriceps dataset only. miR-486 
is encoded within the ankyrin gene and has been shown to 
regulate components of the PTEN/AKT pathway and in doing 
so influences cell cycle kinetics. 40 

In general, microarray and RT-qPCR validation data were in 
close agreement. Those minor differences observed between 
the miRNA microarray and RT-qPCR data are likely to be due 
to the increased sensitivity of the RT-qPCR assays, and to 
the high FDR in the heart array dataset. Methodological dif- 
ferences may be one reason for some of the variations in 
dystromir expression reported previously. 30 ' 3137 We have cal- 
culated PCR efficiencies for all of the TaqMan assays used in 
this study by performing standard curves on serial dilutions 
of cDNA (Supplementary Figure S1a-k and Supplemen- 
tary Table S7). Importantly, in all cases the PCR efficiencies 
were <100%. We therefore utilized the Pfaffl method 41 (which 
corrects for suboptimal PCR efficiencies) to calculate relative 
changes in gene expression as opposed to the Livak (AACt) 
method 42 (which assumes 100% PCR efficiency). Not cor- 
recting for PCR efficiency results in an overestimation of the 
differences between samples and underlines the importance 
of validating assays used in studies of gene expression. 

Multiple lines of evidence point to transcriptional upregu- 
lation of dystromirs in the mdx mouse. Firstly, mature miR- 
NAs that are derived from the same primary transcript are 
concordantly differentially expressed. Secondly, increases in 
the pre-miRNA hairpins are detected in the case of miR-31 
and miR-34c, and lastly, pri-miRNA levels are increased for 
miR-31 , -34c, and -206. A previous report has shown that epi- 
genetic silencing is responsible for changes in the expression 
of certain miRNAs in the mdx mouse as a result of impaired 
nitric oxide signalling and consequent HDAC2 activation. 43 It 
remains to be seen if a similar mechanism is operative for the 
dystromirs we have investigated. 

Several studies have identified miR-1, miR-1 33a and 
miR-206 as being highly enriched (~50-fold) in the serum of 
mdx mice and DMD patients relative to controls. 3233 Intrigu- 
ingly, despite the significant upregulation of these miRNAs 
in serum, this pattern of expression is not reciprocated in 
skeletal muscles (i.e., miR-206 is only increased -4-1 0-fold 
and miR-1 and miR-1 33a are generally decreased or not sig- 
nificantly changed in the skeletal muscles but all three miR- 
NAs are ~50-fold increased in serum). Similarly, the patterns 
of expression of other dystromirs did not generally correlate 
between the muscles and serum. This incongruence between 
the miRNA expression profile of muscle and serum suggests 
that Attdx-enriched extracellular miRNAs do not simply "leak" 
from damaged muscle due to impaired sarcolemmal integrity; 
otherwise the serum profile would be expected to reflect the 
muscle profile. The nonrandom distribution of extracellular 
miRNAs suggests that they may constitute a specific biological 
response. We speculate that miR-1, miR-1 33a, and miR-206 
are actively released from dystrophic muscle and are protected 
from RNase-mediated degradation by either encapsulation in 
microvesicles 44 or in complex with proteins such as Argonate- 
2. 45 We speculate that these extracellular miRNAs function to 
drive muscle regeneration in response to muscle damage or to 
the dystrophic condition through cell-to-cell communication. 
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We have shown that systemic administration of a single 
dose of an antisense oligonucleotide-mediated exon skipping 
therapy resulted in partial normalization of the circulating dys- 
tromirs that are highly enriched in mdx serum. This is consis- 
tent with a previous report that utilized a virus-mediated exon 
skipping strategy. The same study also showed that patients 
with the less severe Becker muscular dystrophy showed inter- 
mediate levels of serum dystromirs between DMD and healthy 
individuals. 32 These serum dystromirs show promise as dis- 
ease biomarkers that could be used to non-invasively monitor 
the effectiveness of experimental therapies in DMD patients. 
Currently, serum creatine kinase (CK) is used as a biomarker 
for muscle damage. In contrast to CK, serum miRNAs have 
been shown to change little in response to exercise 33 and cor- 
relate better with disease severity. 32 A possible limitation of 
this approach is that the levels of serum dystromirs in dystro- 
phic mice and DMD patients are highly variable, as we have 
observed here and in a previous study. 32 This variability may 
mean it is difficult to detect changes in serum miRNA levels 
due to minor improvements in therapeutic outcomes. 

This study has demonstrated muscle-specific changes in 
miRNA expression and indicates that methodological dif- 
ferences are a plausible explanation for the contradictions 
between previous miRNA investigations in the mdx mouse 
and DMD patients. The observation that the miRNA expres- 
sion varies between tissues (even between skeletal muscles) 
is indicative of differential pathological processes occurring 
between these tissues. This unexpected level of complexity 
will need to be considered with regards to the development 
of novel miRNA-based therapies for treating DMD and the 
use of serum biomarkers for monitoring the efficacy of experi- 
mental treatments. 



MATERIALS AND METHODS 

Sample preparation. All animal experiments were carried out 
in the Biomedical Sciences Unit, University of Oxford in accor- 
dance to procedures authorized by the UK Home Office. Male 
C57/B1 0 and mdx mice (n = 4) were sacrificed at eight weeks 
of age. Whole blood was extracted using Microvette CB300 
capillary serum collection tubes (Sarstedt, Leicester, UK) and 
dissection was performed immediately after. The heart, quad- 
riceps femoris, TA, diaphragm and triceps were dissected and 
snap frozen in liquid nitrogen-cooled isopentane. Samples 
were homogenized using Precellys 24 (Bertin Technologies, 
Paris, France) and RNA extracted using TRIzol reagent (Invit- 
rogen, Paisley, UK) as according to manufacturer's instruc- 
tions. Genomic DNA was removed by treatment with TURBO 
DNA- free (Ambion, Warrington, UK) where appropriate. To 
extract total RNA from serum whole blood samples were 
placed at 4 °C for 2 hours and then centrifuged at 13,000 
rpm for 5 minutes in a bench top microcentrifuge. The serum 
supernatant was retained and stored at -80 °C until ready for 
use. Total RNA was extracted from serum using TRIzol LS 
reagent (Invitrogen) as according to manufacturer's instruc- 
tions. RNase-free glycogen (Roche, San Francisco, CA) was 
used as carrier to improve extraction efficiency. 

Pip6e peptide (Ac-RXRRBRRXYRFLIRXRBRXRB-OH, 
with X = aminohexanoyl, B = p-alanine) was synthesized by 



standard Fmoc chemistry and purified by reverse-phase 
high-performance liquid chromatography. The PMO 
(5'-GGCCAAACCTCGGCTTACCTGAAAT) was purchased 
from Gene Tools LLC (Philomath, OR). Pip6e was then conju- 
gated to the PMO through an amide linkage at the 3' end of 
the PMO, followed by purification by cation exchange high- 
performance liquid chromatography, desalted and analyzed by 
MALDI-TOF MS as described elsewhere (C. Betts, A.F Saleh, 
A.A. Arzumanov, S.M. Hammond, C. Godfrey, T. Coursindel et al. 
related manuscript). Pip6e-PMO was dissolved in sterile water 
and filtered through a 0.22-um cellulose acetate membrane 
before use. Pip6e-PMO conjugate was prepared in sterile saline 
for intravenous injection (in a total volume of 160 ul). Under iso- 
flurane anesthesia, 12.5 mg/kg of Pip6e-PMO was administered 
via the tail vein of 12-week-old male mdxm\ce (n = 4). 

miRNA microarray analysis. Implementation and analysis 
of miRNA microarray experiments were designed to comply 
with the minimum information about a microarray experiment 
(MIAME) guidelines. 46 Each RNA sample was analyzed using 
the Agilent 2100 Bioanalyzer (Agilent Technologies, Woking- 
ham, UK) to assess RNA quality. One microgram of total RNA 
was run for each of the 24 samples on Affymetrix GeneChip 
miRNA 2.0 arrays (Affymetrix, Santa Clara, CA) (based on 
miRBase v15). Labeling and hybridization were performed 
according to standard Affymetrix protocols at the Affymetrix 
Core facility BEA, Bioinformatics and Expression Analysis, at 
Novum, Huddinge, Sweden. The processing and data analy- 
ses were performed in Affymetrix miRNA QCTool as well as 
in DNA-Chip Analyzer (dChip) and Excel. The analysis was 
done according to the following workflow; BC-CG background 
adjustment, Quantile normalization and Median Polish summa- 
rization. Since the miRNA 2.0 array covers 131 organisms we 
sorted out Mus musculus transcripts only. Group comparisons 
were performed in dChip software (www.dchip.org). The crite- 
rion for differentially expressed transcripts was set to P < 0.05 
(unpaired, two-tailed Mest). For each comparison, transcripts 
with false detection signals in at least three of four replicates 
were removed. FDR estimations were performed by calculating 
q values (q value R package) from the P value lists. The q value 
of a test measures the proportion of false positives incurred 
when calling that test significant. Array data were curated so 
as to remove ambiguous probe sets and "dead" miRNAs. The 
array data discussed in this publication have been deposited 
in the NCBI Gene Expression Omnibus 47 and are accessible 
through GEO Series accession number GSE36257 (http:// 
www.ncbi.nlm. nih.gov/geo/query/acc.cgi?acc=GSE36257). 

Reverse transcriptase-quantitative PCR. In order to vali- 
date the array results reverse transcriptase-quantitative 
PCR (RT-qPCR) was performed using Small RNA TaqMan 
Assays (Applied Biosystems) as according to manufacturer's 
instructions. All RT-qPCR studies were designed to comply 
with the MIQE guidelines where applicable or practical. 4849 
Briefly, 5 ng of total RNA was reverse transcribed using 
the TaqMan miRNA Reverse Transcription Kit (Applied Bio- 
systems) using miRNA-specific stem-loop RT primers as 
according to manufacturer's instructions. For primary-miRNA 
assays, RNA samples were treated to remove genomic DNA 
using the Ambion TURBO DNase-free kit as according to 
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manufacturer's instructions. 500 ng of DNase-treated total 
RNA were reverse transcribed using the High-Capacity 
cDNA Synthesis Kit (Applied Biosystems) as according to 
manufacturer's instructions. qPCR analysis was performed 
on a StepOne Plus real-time thermocycler (Applied Biosys- 
tems) and TaqMan Gene Expression Mastermix (Applied 
Biosystems). Universal cycling conditions were used (95°C 
for 10 minutes (hotstart) and then 40 cycles of 95 °C for 15 
seconds, 60 °C for 1 minute). For small RNA assays 1.33 ul 
of RT reaction were used in each qPCR. For primary-miRNA 
assays 2 ul of RT reaction were used for each qPCR. All 
primer/probe assays are listed in Supplementary Table S8. 
For mature miRNA assays data were analyzed using the 
Pfaffl method. 41 Gene-of-interest expression was normal- 
ized to miR-16 expression for tissue samples and miR-223 
expression for serum samples as this miRNA has previously 
been used as a reference miRNA in a study of serum dys- 
tromirs in DMD patients. 32 miR-223 expression was highly 
stable across all experimental samples (mean Ct = 25.47, 
SD = 0.28). For primary-miRNA assays data were analyzed 
using relative standard curve method and gene-of-interest 
expression normalized to PpiB expression. 

Reverse transcriptase-PCR. In order to assess the degree 
of exon skipping in treated mdx muscles, 400 ng of total 
RNA was used as a template in a 50 ul RT-PCR using the 
GeneAmp RNA PCR kit (Applied Biosystems). RT-PCR of the 
dystrophin transcript was performed under the following con- 
ditions; 95 °C for 20 seconds, 58 °C for 60 seconds and 72 
°C for 120 seconds for 30 cycles using the following primers: 
DysEx20Fo (5'-C AG AATTCTG CC AATTG CTG AG ) and Dys- 
Ex26Ro (5-TTCTTC AG CTTGTGTC ATCC) . 2 ul of this reac- 
tion was used as a template for nested amplification using 
Amplitaq Gold (Applied Biosystems) under the following con- 
ditions; 95 °C for 20 seconds, 58 °C for 60 seconds, and 72 
°C for 120 seconds for 22 cycles using the following primers: 
DysEx20Fi (5 , -CCCAGTCTACCACCCTATCAGAGC) and 
DysEx26Ri (5'-CCTGCCTTTAAGGCTTCCTT). PCR prod- 
ucts were analyzed on 2% agarose gels. 

Immunohistochemistry. Transverse sections of TA muscle (8 
urn) were collected onto superfrost slides (VWR, Leicester- 
shire, UK). Slides were air dried and soaked in phosphate- 
buffered saline for 10 minutes before being incubated in 20% 
fetal bovine serum and 20% normal goat serum for 2 hours. 
A simultaneous staining protocol was then performed with 
primary antibodies against dystrophin (rabbit polyclonal anti- 
body, Abeam ab15277) and Laminin oc-2 chain (rat mono- 
clonal antibody, Sigma L0663) (20% normal goat serum; 
2 hours) and secondary antibodies Alexa Fluor 488 (goat 
anti-rat IgG) and Alexa Fluor 594 (goat anti-rabbit IgG) (1 
hour). Slides were mounted with Vectashield Hard Set 
mounting medium (with DAPI) (Vector Laboratories, Peter- 
borough, UK). All dilutions and washes were performed with 
phosphate-buffered saline and incubations were performed 
at room temperature. Images were viewed by epifluores- 
cence using a Leica DM IRB microscope and digitally cap- 
tured using Axiovision software (Carl Zeiss Microimaging, 
Jena, Germany). Treated muscle sections were compared 
to samples from untreated age- and sex-matched mdx and 
C57BL10mice (n = 4). 



Statistics. Individual comparisons between disease and wild- 
type samples were tested for statistical significance using a 
two-tailed Student's Mest. For comparisons of more than two 
groups one-way analysis of variance and Bonferroni correc- 
tion post hoc test were performed using GraphPad Prism 5 
(GraphPad Software, La Jolla, CA). Differences were consid- 
ered significant at P values <0.05. 
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Supplementary material 

Figure S1 . Graphs of standard curves for PCR efficiency cal- 
culations for: (a) miR-1, (b) miR-16, (c) miR-21, (d) miR-29c, 
(e) miR-31, (f) miR-34c, (g) miR-1 33a, (h) miR-1 46b, (i) miR- 
199a-3p (j) miR-206, (k) miR-221, and (I) miR-223. 
Table S1. Summary of published DMD-related miRNA profil- 
ing papers. 

Table S2. Table of fold changes, P values and q values for 
differentially expressed miRNAs common to mdx quadriceps, 
diaphragm and heart as measured by miRNA microarray. 
Table S3. Summary of fold changes for all dystromir mdx 
versus wild-type comparisons as measured by small RNA 
TaqMan RT-qPCR. 

Table S4. Validation of miRNA microarray by small RNA Taq- 
Man RT-qPCR. 

Table S5. Concordantly changed mature dystromirs derived 
from common primary-miRNA transcripts in mdx quadriceps, 
diaphragm, and heart. 

Table S6. Concordantly changed mature miRNAs derived 
from common miRNA clusters in mdx quadriceps, diaphragm, 
and heart. 

Table S7. PCR efficiency determinations for RT-qPCR assays 
used in this study. 

Table S8. List of RT-qPCR assays used in this study. 
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